Current-induced noise and damping in non-uniform ferromagnets 
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In the presence of spatial variation of the magnetization direction, electric current noise causes a 
fluctuating spin-transfer torque that increases the fluctuations of the ferromagnetic order parameter. 
By the fluctuation-dissipation theorem, the equilibrium fluctuations are related to the magnetiza- 
tion damping, which in non-uniform ferromagnets acquires a nonlocal tensor structure. In biased 
ferromagnets, shot noise can become the dominant contribution to the magnetization noise at low 
temperatures. Considering spin spirals as a simple example, we show that the current-induced noise 
and damping is significant. 

PACS numbers: 72.70.+m, 72.25.Mk, 75.75.+a 
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Electric currents induce magnetization dynamics in 
ferromagnets. Three decades ago, Berger— » showed that 
an electric current passing through a ferromagnetic do- 
main wall exerts a torque on the wall. The cause of this 
spin-transfer torque is the reorientation of spin angular 
momentum experienced by the electrons as they adapt to 
the continually changing magnetization. Subsequently, it 
was realized that the same effect may also be present in 
magnetic multilayers 3 -. In the latter case, the torque may 
cause reversal of one of the layers, while in the former, 
it may cause domain wall motion. The early ideas have 
been confirmed both theoretically and experimentally 4 . 

Recently, the importance of noise in current-induced 
magnetization dynamics has drawn attention. Although 
often noise is undesired, it may in some cases be quite 
useful. Wetzels et al& showed that current-induced mag- 
netization reversal of spin valves is substantially sped up 
by an increased level of current noise. The noisy cur- 
rent exerts a fluctuating torque on the magnetization 6 . 
Ravelosona et al. 1 reported observation of thermally as- 
sisted depinning of a narrow domain wall under a cur- 
rent. Thermally- assisted current-driven domain wall mo- 
tion has also been studied theoretically^. 

The present paper addresses current-induced magneti- 
zation noise in non-uniformly magnetized ferromagnets. 
The spatial variation of the magnetization direction gives 
rise to increased magnetization noise; by a fluctuating 
spin-transfer torque, electric current noise causes fluc- 
tuations of the magnetic order parameter. We take 
into account both thermal current noise and shot noise, 
and show that the resulting magnetization noise is well 
represented by introducing fictitious stochastic magnetic 
fields. By the fluctuation-dissipation theorem (FDT), 
the thermal stochastic field is related to the dissipation 
of energy, or damping, of the magnetization. The FDT 
hence constitutes a simple and efficient way to evaluate 
the damping, providing also a physical explanation in 
terms of current noise and spin-transfer torque. Since the 
correlator of the stochastic field in general is inhomoge- 
nous and anisotropic, the damping is a nonlocal tensor. 
As a simple and illuminating example we consider ferro- 



magnetic spin spirals, for which the field correlator and 
damping become spatially independent. It is shown that 
for spirals with relatively short wavelength (~ 20nm), 
the current-induced noise and damping is substantial. 
Since half a wavelength of a spin spiral can be consid- 
ered as a simple model for a domain wall, this suggests 
that current-induced magnetization noise and damping 
should be an issue for narrow domain walls. 

It is instructive to start with an introduction to the 
FDT for uniform (single-domain) ferromagnetic systems, 
characterized by a time-dependent unit magnetization 
vector m(t) and magnetization magnitude M s (the sat- 
uration magnetization). The spontaneous equilibrium 
noise of such macrospins is conveniently described by the 
correlator Sij(t — t') = (5nii(t)8mj(t')) , where 5nii(t) — 
m i{t) — (^j(^)) is the random deviation of the magne- 
tization from the mean value at time t. The brackets 
denote statistical averaging at equilibrium, while i and j 
denote Cartesian components. The magnetization fluc- 
tuations are assumed weak, so that they to first order are 
purely transverse to the equilibrium (average) direction 
of magnetization. Applying a weak external magnetic 
field h(° xt )(t), the magnetization can be excited from the 
equilibrium state. Assuming linear response, the result- 
ing transverse change in magnetization is 



Ami(t) 
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dt'xijit-t'^it'), 



(1) 



defining the transverse magnetic susceptibility XijiP ~ 
as the causal response function. The FDT relates this 
susceptibility to the equilibrium noise correlator—: 



Sn (t-t') = 



M S V 
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where T is the temperature and V is the volume of the 
ferromagnet. 

The spontaneous equilibrium fluctuations 5m(t) may 
be regarded to be caused by a fictitious random magnetic 
field h(t) with zero mean. We can derive an alternative 
form of the FDT in terms of the correlator (hi(t)hj(t')). 
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To do so, simply note that Eq. JTJ) implies that <Jm,(u;) = 
yV Xij(u))hj(u)) in Fourier space. Inverting this relation, 
it follows from Eq. ^) that 



<MW)) 



xt,V) 



(3) 

where Xy ( w ) is the ij-component of the Fourier trans- 
formed inverse susceptibility tensor. 

The magnetic susceptibility can be found from the 
Landau-Lifshitz-Gilbert (LLG) equation of motion. The 
stochastic LLG equation describes magnetization dynam- 
ics and noise in both uniform as well as non-uniform fer- 
romagnets, and reads 

^ = - 7 mx[H off + h + h( cxt )]+a mx^. (4) 

Here 7 is the gyromagnetic ratio, H c ff is an effective static 
magnetic field determining the equilibrium state, h(t) is 
the above random noise- field, h( ext ) (t) is the weak excita- 
tion introduced in Eq. {T]), and ao is the Gilbert damping 
constant. Linearizing this equation in the magnetic re- 
sponse to h( ext ) (t) , we find the inverse susceptibility 
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7|H off | - iuia Q iu 

~ilu 7|H cff | - iujao 



(5) 



written in matrix (tensor) form in the plane normal to the 
equilibrium magnetization direction. Note that the static 
field has here been assumed local and magnetization in- 
dependent. While not valid in most realistic situations, 
this simple form for the effective field captures the key 
physics of interest here, since only the dissipative part 
of the susceptibility (the Gilbert damping term) affects 
the noise. Inserting Eq. © into Eq. ([3]), we get the 
well-known result^ 
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where i and j denote components orthogonal to the equi- 
librium magnetization direction. This expression relates 
the equilibrium noise, in terms of h, to the damping or 
dissipation of energy in the ferromagnet. It may be noted 
that in thin ferromagnetic films in good electrical contact 
with a metal, the equilibrium noise and corresponding 
Gilbert damping has been shown to be substantially en- 
hanced. This is due to the transfer of transverse spin 
current fluctuations in the neighbouring metal to the 
magnetization-^. 

We now turn our attention to a more complex sys- 
tem, i.e., a metallic ferromagnet whose direction of mag- 
netization m is varying along some direction in space, 
say, the y-axis. It is assumed that the spatial variation 
is adiabatic, i.e., slow on the scale of the ferromagnetic 
coherence length. The ferromagnet is furthermore as- 
sumed to be translationally invariant in the x- and z- 
directions, and its magnetization magnitude is taken to 
be constant and equal to the saturation magnetization 



M s . In general, the dynamics and fluctuations of such 
a magnetization texture depend on position. Due to the 
spatial variation of the magnetization, longitudinal (i.e., 
polarized parallel with the magnetization) spin current 
fluctuations transfer spin angular momentum to the fer- 
romagnet. The resulting enhancement of the magnetiza- 
tion noise is described by introducing a random magnetic 
field, whose correlator is inhomogenous and anisotropic, 
unlike Eq. ([6]). By the FDT, the correlator is related 
to the magnetization damping, that acquires a nonlo- 
cal tensor structure. In the following we make use of 
the fact that the time scale of electronic motion is much 
shorter than the typical precession period of magneti- 
zation dynamics, as implicitly done already in Eq. |6]). 
We shall disregard spin-flip processes and the associated 
noise. Spin-flip corrections in Fe, Ni, and Co are expected 
to be small because the spin-flip lengths are long com- 
pared to the length scale of spatial variation (domain 
wall width) we consider. Spin-flip is important in Py. 
However, domain walls in Py are so wide that the effects 
discussed here are not important anyway. We therefore 
do not discuss spin-flip scattering. 

It is convenient to transform the magnetization tex- 
ture to a rotated reference frame, defined in terms of the 
equilibrium (average) magnetization direction mo(j/) = 
(m(?/,i)) of the texture. The three orthonormal unit 
vectors spanning this position-dependent frame is vi = 
v 2 x v 3 , v 2 = (dm /dy)/\dm /dy\ and v 3 = m . The 
local gauge 



U{y) = [ vi(y) v 2 (y) v 3 (y) 



(7) 



transforms the magnetization, and hence the relevant 
equations involving the magnetization, to this frame. 
That is, Umo(y) = mo = z, where the tilde indicates 
a vector in the transformed frame. We note also that 
Uvi = x and UV2 — y, and that U is orthogonal, i.e., 

U- 1 =U T = [V! v 2 v 3 ]. 

We consider a charge current I flowing through the 
ferromagnet along the y-axis. Assuming that the equi- 
librium magnetization direction mo(y) changes adiabat- 
ically, the electrons spins align with the changing mag- 
netization direction when propagating through the tex- 
ture. The spin current is then anywhere longitudinal, and 
hence given by l s (y) — I s mo(j/). The alignment of the 
electrons spins causes a torque t(jj) = dl s (y)/dy on the 
ferromagnet. Since dl s (y)/dy is perpendicular to mo(y), 
the torque can be written r(y) = — mo(j/) x [mo(y) x 
dl s (y)/dy], orf(y) = Ur(y) = -m x [m x Udl s (y)/dy] 
in the transformed representation. When 1 = 0, which 
we will take in the following, /„ = and f = 0, on aver- 
age. However, at T ^ thermal fluctuations of the spin 
current result in a fluctuating spin-transfer torque 

AT(y,t) = -A/,(t)mox[moxl7^M] ( 8 ) 

dy 

where A/ s (t) are the time-dependent spin current fluctu- 
ations with zero mean, propagating along the y-direction. 
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The action of the fluctuating torque on the magnetiza- 
tion is described by the LLG equation if we, by conserva- 
tion of angular momentum, add the term ^At/^MsA) 
on the right hand side. Here A is the cross section 
(in the xz-plane) of the ferromagnet. Linearizing and 
transforming the LLG equation to the rotated reference 
frame, it is seen that the fluctuating torque JHJ can 
be represented by a random magnetic field h'(y, t) = 
AI S (t)/M s A) [mo x Udmo(y)/dy], analogous to h(t) dis- 
cussed above. Using Eq. ([7]) 



AI s (t) 
M s A 



dm (y) 



dy 



(9) 



i.e., the (transformed) current-induced random field 
points in the ir-direction. 

The longitudinal spin current fluctuations AI s (t) can 
be found by Landauer-Biittiker scattering theoryS^. 
Disregarding spin-flip processes, the spin-up and spin- 
down electrons flow in different and independent chan- 
nels. In the low- frequency regime, in which charge is in- 
stantly conserved, longitudinal spin current fluctuations 
are perfectly correlated throughout the entire ferromag- 
net. Hence, the thermal spin current fluctuations are 
given by&i^ 

(AI s (t)AI s (t>)) = ^2fc B T(G T + GMt t>), (10) 

where Gfn) is the conductance for electrons with the spin 
aligned (anti)parallel with the magnetization. This ex- 
pression is simply the Johnson-Nyquist noise generalized 
to spin currents^. We find from Eqs. ((9|) and (fPTJf 



(h' x (y,t)K(y',t>)) = 2kBT ^ y,) 5(t - (ii) 



for the correlator of the current induced random field. 
Here we have defined 
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dm (y) 



dy 



dm (y') 



dy 



(12) 



with a = (Gf + G{)L/A the total conductivity. Recall 
that h' y (t) = h' z (t) = 0. Eq. (fl~Tj) describes the nonlo- 
cal anisotropic magnetization noise due to thermal cur- 
rent fluctuations in adiabatic non-uniform ferromagnets. 
This excess noise vanishes with the spatial variation of 
the magnetization. As a consequence of Eq. ([lT)|) , the 
random field correlator depends nonlocally on the mag- 
netization gradient. 

According to the FDT, the thermal noise is related to 
the magnetization damping. Since the noise correlator 
(llip is inhomogeneous and anisotropic, the correspond- 
ing damping must in general be a nonlocal tensor. To 
evaluate the damping, we hence need the spatially re- 
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FIG. 1: An example of a non-uniform ferromagnet. The mag- 
netization rotates with wavelength A in the ?/2-plane, forming 
a spin spiral. 



solved version of the FDT, which reads 
(6m i (y,t)8rh j (y' ,t')) 



M s A 

Xij(y,y',v) - Xji(y',y,y) 



i2nu) 



(13) 



in the transformed representation. Here 8rh(y, t) — 
USm(y,t) — 8m x (y,t)x + 8m y (y,t)y are the spatially 
dependent transformed magnetization fluctuations. The 
susceptibility is defined by 



Arhi{y,t) = 



j 



dy'dt'xi^y'^-t'fh^iy'J), 



analogous to Eq. fl} 
and magnetic excitations transformed: h 



(14) 

but with the external field 
~^\y,t) = 



Uhf**\y,t) and Am(y,t) = UAm(y,t). The suscep- 
tibility in the local gauge frame differs from Eq. (J5J) and 
has to be determined. It is straightforward to gener- 
alize Eqs. (fTB"]) and (|14p to the case of general three- 
dimensional dynamics. 

We may substitute hf xt) (y', t') by h'j (y' , t') in Eq. JHJ 
to find the fluctuations 8m(y, t) of the magnetization vec- 
tor caused by the spin-transfer torque. Combining this 
expression with Eqs. (|13|) and pip , we arrive at an inte- 
gral equation for the unknown susceptibility, from which 
the nonlocal tensor damping follows. Instead of finding 
a numerical solution for an arbitrary texture, we con- 
sider here a ferromagnetic spin spiral as shown in Fig. [TJ 
for which the description of magnetization noise can be 
mapped onto the macrospin problem. A simple analytical 
result can then be found, allowing for a comparison with 
Eq. (O, and hence an estimate of the relative strength 
and importance of the current-induced noise and damp- 
ing. 

Spin spirals can be found in some rare earth metals^ 
and in the 7-phase of ironi^, and are described by 
m o(y) = [0,sin%),cos%)], where %) = 2?ry/A = qy, 
with A the wavelength of the spiral. Then dm^(y) / dy = 
q[0, cos#(y), — sm9(y)] so that \dm (y)/dy\ = q. As em- 
phasized earlier, our theory is applicable when the wave- 
length is much larger than the magnetic coherence length. 
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For transition metal ferromagnets, the coherence length 
is of the order of a few angstrom. From Eq. (TT2"|) we 
find £ xx = ^T^aq 2 /(4e 2 M s ). The current-induced noise 
correlator (jlip for spin spirals is hence homogeneous, 



jM s V 



-5(t-t% 



(15) 



similar to Eq. ©, but anisotropic. The problem of relat- 
ing noise to damping in terms of the FDT can therefore 
be mapped exactly onto the macrospin problem: The 
transformation Q can be used to show that equations 
analogous to Eqs. CP)-© are valid for the spin spiral, 
when analyzed in the local gauge frame. It is then seen 
that the damping term corresponding to Eq. (| 1 5j) is 



in the transformed representation. Here 







(16) 



(17) 



is the 2x2 tensor Gilbert damping in the xy-plane. 
Hence, £ xx is the enhancement of the Gilbert damping 
caused by the spatial variation of the magnetization and 
the spin-transfer torque. Due to its anisotropic nature, 
£ is inside the cross product in Eq. (|16p . ensuring that 
the LLG equation preserves the length of the unit mag- 
netization vector rh. 

In order to get a feeling for the significance of the 
current-induced noise and damping, we evaluate £ nu- 
merically for a spin spiral with wavelength 20 nm, and 
compare with cto. Taking parameter values for cto, M s , 
and a from Refsi 16 ' 17 i 18 ' 19 , we find £ xx « 5ao for Fe (with 
a = 0.002), and £ xx « 4a for Co (with a Q = 0.005). 
Hence, current-induced noise and damping in spin spi- 
rals can be substantial. Considering half a wavelength 
of the spin spiral as a simple domain wall profile, these 
results furthermore suggest that a significant current- 
induced magnetization noise and damping should be ex- 
pected in narrow (width ~ 10 nm) domain walls in typ- 
ical transition metal ferromagnets. The increased noise 
level should assist both field- and current-induced do- 
main wall depinning£i2i2£. The increased damping should 



be important for the velocity of current-driven walls, 
which recent theoretical and experimental advances sug- 
gest is inversely proportional to the damping 4 -. The in- 
creased noise and the tensor nature of the Gilbert damp- 
ing should be taken into account in micromagnetic sim- 
ulations. 

So far we have only considered thermal current noise; 
let us finally turn to shot noise. With the voltage U 
across the ferromagnet turned on, a nonzero current / 
flows in the y-direction. Disregarding spin-flip processes, 
the resulting spin current shot noise i a 6 ' 13 



(AI^\t)AI^\t>)) 



(2e) 



-eUFG5(t - 1') 



(18) 



at zero temperature. Here the superscript (sh) empha- 
sizes that we are now looking at shot noise. The Fano fac- 
tor F is between and 1 for non- interacting electrons- 21 . 
When the length of the metal exceeds the electron- 
phonon scattering length X ep , shot noise vanishe d 13 ' 21 . 
A ep is strongly temperature dependent, and can at low 
temperatures exceed one micron in metals. To find the 
contribution from shot noise to the magnetization noise, 
simply replace Eq. (fTT)|) with Eq. (TT5| in the above cal- 
culation of the random-field correlator. In experiments 
on current-induced domain wall motion, typical applied 
current densities are about j = 10 s A/cm 2 —. At low 
temperatures, the ratio of shot noise to thermal current 
noise, eUF/2kBT, can then exceed unity for long (but 
not longer than X ep ) ferromagnetic (e.g. Fe) wires. Shot 
noise can hence be expected to be the dominant contri- 
bution to the magnetization noise at low temperatures. 

In summary, we have calculated current-induced mag- 
netization noise and damping in non-uniform ferromag- 
nets. Taking into account both thermal and shot noise, 
we evaluated the fluctuating spin-transfer torque on the 
magnetization. The resulting magnetization noise was 
calculated in terms of a random magnetic field. Em- 
ploying the FDT, the corresponding enhancement of the 
Gilbert damping was identified for spin spirals. 
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